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Abstract 

A formally exact set of equations is derived for the description of nonequi- 
librium phenomena in classical liquids and glasses. With the help of a non 
equilibrium projection operator formalism, the correlation functions and fluc- 
tuation propagators are expressed in terms of memory functions and time 
dependent collective frequencies. This formally exact set of equations is ap- 
proximated by applying mode coupling approximations to the memory func- 
tions. The resulting set of equations for wavevector dependent correlation 
functions, fluctuation propagators and one-time structure factors Sq(t) gener- 
alizes the well known mode coupling theory of the glass transition to situations 
far away from equilibrium. 
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I. INTRODUCTION 



The mode coupling theory of the glass transition (MCT) is based on the experimental 
and theoretical understanding of dense liquids. Although approximative in nature, mode 
coupling theories were quite successful to describe quantitatively the dynamics of dense 
liquids close to the triple point. The dynamics of dense liquids is dominated by two mecha- 
nisms, the precursor of the cage effect and a complicated coupling of current fluctuations and 
density fluctuations. This coupling is the main cause of activated motion close to the glass 
transition [|T],0 and will not be considered further in this paper. The cage effect is a collective 
phenomenon caused by the mutual hindering of the motion of the atoms or molecules due 
to the strong repulsion of the cores of the atoms. Instead of a unidirected flow-like motion 
a kind of stop and go mechanism is observed in simulations of dense liquids. The particles 
tend to rattle in a cage formed by the surrounding atoms. To escape this cage an atom has 
to find a hole in the surrounding wall formed by the other atoms. The same is of course true 
for the atoms building the wall. Due to the thermic motion a situation will arise in which the 
particles can pass each other. This process will be repeated over and over again. In this way 
the motion of particles, which are nearby at some instant of time, will decorrelate slowly. In 
the process a back flow pattern is created. By increasing the density, the local cage effect 
will eventually lead to a complete blocking of long range motion of the particles. They will 
be stuck in their local cages and thus form a glass. The details how this is done, especially 
how this glassification affects the various spatial scales are described by the idealized mode 
coupling theory Within this theory it is e.g. possible to understand the effect of the 
stiffening of the local cages on the thermodynamics and hydrodynamics close to the glass 
transition [Q. Experimentally observed phenomena like the frequency dependent specific 
heat are thus explained quite naturally in a unified framework. 

The underlying mathematical mechanism of the glass transition in the idealized MCT is 
a bifurcation in the equation for the dynamics of the density - density correlation function 
of a system of N particles. Using a standard Mori projection operator formalism for the 
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density and current fluctuation nq{t),jq{t) with 



N 



n, 



,{t) = exp{iLt)nq{0) = ^exp{iqfi{t)) (1) 



- N ^ 

jq{t) = -jq{t) = Yl -Viexp{iqfi{t)), (2) 

Q i Q 

the following formally exact equation for the density density correlation function can 
be derived: 

d(j)(t) , , . , 

= ^(iKl^t) ■= Mn,3 (3) 

Here the correlation function is given as 

0(t) = (/)„„(t) = {nq\nq{t)) := Tr{exp{-j3H)n*q{Q)nq{t)) (5) 

L is the Liouville operator. Its action on some dynamic variable is defined by the Poisson 
bracket 

N 

LA = -iY,{H,A}. (6) 

i 

By construction the memory function M is free of any hydrodynamic singularity and cannot 
be expressed as a linear functional of the correlation functions. Instead it can be written as 
a correlation function of fluctuation forces. 

Mqit) = {QLjql exp{iLQt)QLjq)/S,,. (7) 

Sjj is the static current current fluctuation function. In equilibrium it is given as 

Sjj = NkBT/m (8) 

where m is the mass of an atom and T is the equilibrium temperature. In the mode cou- 
pling approximation the fluctuating forces are approximated as a superposition of two mode 
contributions [HI 



QLjg ^^v{q\kp)nknp (9) 

kp 

v{q\kp) = iQLjq\nknp){SkSpy^ = 5^^^_j:Sjjv{k, q - k) (10) 

The second equality in eq. (0) is due to momentum conservation. Sq is the static structure 
factor (j){q,t = 0). The vertex v can be evaluated in equilibrium. Neglecting static direct 
three point correlation functions the memory function is given by |^ 

Mq{t) = n{kBT/mf\ ^'i^, q - k)<P{q - k, t)(P{k, t) (11) 
^ k 

with 

— # — * 

vik,q- k) = {khk + {q- k)^Cq.ky (12) 

Here Cg is the (equilibrium) direct correlation function. Eqn |^ and |Tl| form a closed set 
of equations for the density - density correlation function in equilibrium. The solution 
properties of the above equilibrium MCT (EMCT) equation are very well known They 
exhibit an ideal glass transition as a bifurcation scenario from an ergodic to a non ergodic 
state. The transition temperature Tc (or density, pressure ...) is independent of the cooling 
rate since it is an equilibrium property of the system. By approaching the structural 
relaxation is dramatically slowing down, developing the typical two step relaxation pattern 
of /3 - and a - relaxation. The EMCT was quantitatively tested for a experimental realization 
of a hard sphere system . Without any fit factor (up to an overall scale factor) it is possible 
to reproduce with 20 % accuracy the density - density correlation for the wavevector at the 
peak of the structure factor. With the same kind of approximations for memory functions 
the theory was also quantitatively tested for a binary Lennard Jones systems [Q, a system 
of dumbbells and a model system for water Any generalization to non equilibrium 



situations should therefore in equilibrium reduce to the above set of equations (|) and (p^. 

The non equilibrium situations to be described in this paper arise if one investigates 
the relaxation in a system, which is not yet equilibrated. When cooling down from some 
temperature Ti > to T2 > T^, one has to wait at least for some time comparable 
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to a typical a - relaxation time TqiT'i) to guarantee equilibration at the temperature T^. 
As indicated, the a - relaxation time depends on the wavevector q. To use the EMCT, 
density fluctuations for all wavevectors have to be equilibrated. For hard sphere like systems 
the largest a - relaxation time is given by the wavevector at the peak of the structure 
factor i.e only for a time tyj > r^qmax) the above equations are able to describe the density 
correlation functions. For times smaller than T{qmax) nonequilibrium effects like breaking of 
time translation invariance are to be expected. E.g. the one-time structurefactor Sq{t) = 
(pq(t,t) is time dependent, since the probability density has not yet reached its equilibrium 
value. 

As soon as the system is quenched below T^, the behavior of the system changes qual- 
itatively. In EMCT the a relaxation time is going to inflnity for all wavevectors, when T 
is approaching from above, i.e. the system never reaches equilibrium for T < T^. This 
means that e.g. even for inflnite long waiting time the one-time structure factor will be dif- 
ferent from the value, it would reach in the presence of hopping processes. For real systems, 
where activated processes cannot be neglected, the waiting time has to be smaller than the 
inverse of the typical hopping rate, in order to see aging. Therefore the applicability of 
EMCT has to be reconsidered within the larger framework of a theory, which is applicable 
also in nonequilibrium situations. 

In the theory of spin glasses quite some progress in this direction was obtained in recent 



years. Already Kirkpatrick, Thirumalai |jTT| showed that the dynamics of an infinite ranged 
p - spin system with p > 2 and a disordered Potts system {p > 4) exhibit great similarities 
to the dynamics of structural glasses at least close to the critical temperature (which 
is in the spin glass literature called To). The equation for the correlation function for 
T > Tc reduce to a schematic EMCT model of the glass transition. For p = 3 the p - spin 
model is equivalent to the so called Leutheusser model [0] above T^. This model exhibits a 
glass transition, but its a - relaxation is a simple exponential relaxation and can therefore 
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not explain the typical stretched relaxation observed in real glasses. Q In the context of 
superccoled liquids the Leutheusser model can only be motivated by the complete neglect 
of the wave vector dependence in equ. For the p - spin model the Leutheusser model 

is an exact theory due to the mean field character of this system. Cugliandolo and Kurchan 
showed for the p - spin model , that the EMCT is not valid any more below Tc, due to 



similar reasons as explained above [|1^]. Instead they derived nonequilibrium mode coupling 
equations for this model and solved them asymptotically below Tc. Contrary to the EMCT, 
the dynamics cannot be understood in terms of correlation function alone, but a coupled set 
of equations for the correlation functions 0(t^,t) and susceptibilities is necessary. 

The most interesting property of the solution of these equations is the separation of the 
dynamics in a short time part, where the fluctuation dissipation theorem (FDT) is obeyed 
and an aging part where FDT is violated in a very specific way. In the aging regime the 
equations for t), t) can be reduced to a single equation for a function 0(A) with 
A = h{tu,)/h{t). The susceptibility is given as 

xK.t) = ^m (13) 

„ , , , X d(j) 
^^^^ = TdX 

The function h has to be determined by matching the FDT solution and the aging regime. 



So far its form is only known for the spin glass with p = 2 [jl4[. The term - was interpreted 



as an effective temperature [15|. Since the important work of Cugliandolo and Kurchan a 



large number of papers has addressed the problem of nonequilibrium relaxations in mean 



field models |T6|. It was shown for general nonlinear fieldtheories that in the cases where 
NMCT equations are exact, a mapping to an equivalent disordered problem can be found 
||17| |. This can also be achieved by deterministic but highly irregular interactions [jTB . 



^It is possible to consider a more general class of schematic models, which show non exponential 
relaxation B, but they are not easy to justify within the theory of liquids. 
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Simulation of binary liquids in the out of equilibrium regime |l9|-[2^ exhibit striking 
qualitative similarities to the dynamics found in the p -spin model. There is the separation 
in FDT and aging regime. Also the scaling behaviour (|1^) seems to be fulfilled. It is 
worth mentioning, that the crossover from the FDT regime to the aging regime, i.e. the 
behaviour of the function x{tw,t), seems to be different than in the p -spin model. To 
be able to describe behaviour of real liquids the schematic models are not sufficient, since 
they lack the important wavevector dependence, leading e.g. to non exponential relaxation 
and experimentally testable results like wavevector dependent Kohlrauschexponents |^ 



Equations for multicomponent models have been formulated on an abstract level fT^,^, 
but there relation to the theory of glass forming liquids is not obvious. Additionally there 
are no equations for the one time structure factor Sq{t) for real liquids. In the following, I 
am going to formulate a generalization of the EMCT for real liquids, which will enable the 
study of the dynamics of real liquids above and below Tc in equilibrium and nonequilibrium 
situations. 



II. THE NONEQUILIBRIUM PROJECTION OPERATOR FORMALISM 

The main reason for the appearance of generic nonequilibrium relaxation is the inability 
of the system to equilibrate below (if hopping is neglected, which I will assume in the 
following). To keep the formalism as simple as possible, I will consider a situation where the 
probability density pNE is different from the canonical (p 7^ exp(— or other equilibrium 
densities but the microscopic dynamics is still Hamiltonian. The only requirement for not 
being in equilibrium is 

{PNE, H} ^ (14) 

This situation may e.g. arise if the system is isolated after an initial quench from an 
equilibrated liquid above Tc to a glass below Tc (e.g. by coupling to an external heat bath, 
whose temperature is below ). The system tries to get into equilibrium with the heat bath. 
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As argued in it can be expected that the kinetic energy is equihbrating on microscopic 
time scales. It is the potential energy which does not reach its equilibrium value due to its 
contribution from structural relaxations. After the initial quench we may therefore expect 
to describe the system with an initial density operator pne = exp(— (/3/i^' + /5/e^°*). Here /?/ 
is the inverse of the initial temperature before the quench and /5/ is the inverse of the final 
temperature. For Pc = ^/{kBTc) the relation [3i < j3c < Pj is fulfilled, if, e^°* is the kinetic 
and potential energy respectively. If needed the formalism can also be developed by assuming 
a coupling to an external heat bath by using a Gaussian thermostat to guarantee the average 
of the kinetic energy being constant at any time ||2^, but here I want to concentrate on the 
more simple case of an isolated system described by a nonequilibrium density operator. 

As basic variables for the following projection operator formalism (POF) Ai{t) I again 
will use density fluctuations (P and current fluctuations (0). The energy fluctuations are 
neglected. They are of course influenced by the glass transition, but they do not qualitatively 
change the glass transition scenario |Q. The main idea of any projection operator formalism 
is to separate the dynamics of the variables in a part directly proportional to a chosen 
variable, a part which is given by the fluctuation of the variables in the past and a fluctuating 
force, whose dynamic was never directly proportional to the fluctuations of the basic variable 
6Ai{t). The most direct generalization of the standard Mori formalism to situations, in which 



also equal time correlation functions are time dependent is obtained by using p7H29|| 



P{t) = exp{-iLt)6Ai{0)S-^\t)6Ak{t)exp{iLt) (15) 

:= exp{-tLt)\Ai{0))Sr,\t){Akit)\exp{tLt) 
Q{t) = 1 - Pit) (16) 

Equal indices are summed over. In equilibrium this projection operators reduces to the 
time independent Mori projection operator in equilibrium. The matrix S is the equal time 
correlation matrix of the variables Ai 

Sik = {Ai{t)\Ak{t)) = Tr{pME exp(zLt)A,(0)| exp(2Lt)Afc(0)). (17) 



The projection operators have the following properties. 

P{ti)P{h) = P{h) 

Q{ti)Q{t2) = Qih) (18) 
Q{t2)P{h) = 0. 

With the help of (|15|,|T6|) a formally exact set of equations for the fluctuation of the 
density (|I|) and the current (|^) can be derived 

dnq{t) 



dt 
djq{t) 



-- ^q3,{t) (19) 

dt '^^^ = (2°) 

i5nq{t)q^Vtnj{t) + i5jq{t)q^VLjj{t) 

- f dt'6nq{t')S-}{t')Mq{t',t) - f dt'6jq{t')Sjl{t')Mq{t',t) 

+i exp{iLtj^)Q{t^)G{t^, t)L6jq{0) 

The structure of the equation is very similar to the equilibrium equation. Eq. (|1^) is 
simply the continuity equation for the density. q^Vtnjj'f^jj are generalized time dependent 
frequencies. They are given by 

^n, = S-l(t){nq(t)\L3q{t))/q + S-lmq(t)\Ljq[t))/q (21) 
= Sj^{t){nq{t)\Ljq{t))/q + Sj^{t){jq{t)\Ljq{t))/q (22) 

In equilibrium Vtjj would be zero due to time translational invariance, but has to be kept 
in nonequilibrium for finite times. The memory function Mq is an explicit function of two 
times. 

Mq{ti,t) = iLjq{h)\exp{iLt^)Q{t^)Git^,t)L5jq{0)) (23) 
G{ti,t) = Texpit t dt' LQ{t')) (24) 

Where the operator T induces a time ordering from left to right. It is important to note 
that in general the memory function (|23|) does not have the form of a correlation function 
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of fluctuating forces, but is the correlation of the accumulated random force between ti and 
t and the real force Ljq{ti) at time ti. It describes the dissipation between ti and t ^8 



This will be important for finding appropriate approximations. The last term in (|20| ) is a 
fluctuating force Fqit^.t) 

Fg(t^, t) = exp{iLtMtn^)G{tvj, t)Ljg{0). (25) 



Due to the properties of the projection operators ([T8|) the time ordered product in (^31) 
propagates the force Ljq always in the space of functions perpendicular to density and current 
at every time step. The fluctuating force Fq{tw, t) is not correlated with the fluctuation of 
density or current at time t^. 

Equ. (|19|, |20|) have the form of an inhomogeneous integro differential equation. They can 
be solved formally by introducing fluctuation propagators N^i, which obey the homogeneous 
part of equ. (H, l^). 

dNin{q;ti,t) 

— = tqNij (g; 1 1 , t ) (26) 



dNij{q;ti,t) 



iNi^{q- ti, )^]„,(t) + iNijiq; t,,t)njj{t) (27) 
dt'N,r,{q;h,t')S-}{t')Mq{t',t) 
dt'N,,{q;h,t')Sjl{t')Mq{t',t) 



dt 

ft 



The initial conditions are Nik{q]t,t) = 6ik- If a density or current fluctuation is given at a 
time ti e.g. by applying a known external field to the system, the fluctuations at all later 
times are then given given by 

6nq{t) = 6Ai{ti)Ni4q; h, t) + i f Fg{ti,t')N,n{t' , t) (28) 
Sjq{t) = 6Ai{ti)Ni,{q; h, t) + i f F,(ti, t')N,j{t\ t), (29) 



where Fg{ti,t) is the random force eq. (^5]). The Ai are density and current respectively for 
i = 1,2. This can easily be checked by replacing density and current in equ. ( p!9| , pOD by (|28| , 
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The fluctuation propagators are very closely related to the susceptibilities. These are 
jiven by Poisson brackets of the variables |^ . 



Xik = {{SA:ih),5A,m (30) 
The density susceptibihty can therefore be written with the help of eq. ( P^ ) as 



Xnn(tl, t) = -tqNjn{tl,t) + % / ({^^(ti), F,(ti , t') } A^in(t', t)) (31) 

Here it was used that Xnra(^i;'^i) = and Xnj{(l\t\,t\) = —iq. In mode coupling approxi- 
mation, the second term in (^) is zero The dynamical susceptibility can therefore be 
identified with the fluctuation propagator —iqNjn{q;ti,t). 

The equation for the density correlation function can be easily derived from eqs. 
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— -^^ — = q ^nj[t) + 1 — — qS^jA^) 

-iq dt'(l),{h,t')S-}it')M,{t', t) (32) 
+ f\-iqNjn{q]t' M))M,{t' ,t) 







The function Mq{t\ t) is in general not identical to the memory function Mqiti^t) in (p3[), 
but it is the correlation function of the fluctuating forces. 

Mq{h,t) = (g(o)F,(o,ti)|g(o)F,(o,t)) (33) 



Eqs. (|3^j27|) are a formally exact set of equations for density correlation function and 
fluctuation propagators, which reduces to the response or integrated response within mode 
coupling theory. If we assume equilibrium, i.e. if we use the canonical ensemble, the set 



of equations reduces exactly to the one in equilibrium p6|. Time translation invariance 



is restored, the fluctuation propagator Nnn(ti,t) is identical to the normalized correlation 
function, M = M and the equation (|32D reduces to (H) with q'^flnj = 
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For nonequilibrium situations a further step has to be taken. In EMCT the only input 
needed is the static structure factor. The frequency Qg can easily be calculated if S{q) is 
known. In nonequilibrium S{q) turns into a time dependent function. E.g. for quenches we 
only know the value of Sq{t = 0), which is just the equilibrium structure factor of the liquid 
at the begin of the quench. The equations presented so far, do not provide an equation for the 
equal time correlation function. It is easily possible to obtain an equation by differentiating 
twice Sq(t). 

fj'^Q (f) 

= 2q\jMm) - 2Re{{nMLj{t)) (34) 

Equations like this cannot be treated with the standard non equilibrium projection op- 
erator formalism, since it is assuming implicitly that frequencies like q{nq\Ljq) are known 
or can be simply calculated. This is not the case for liquids close to the glass transition in 
nonequilibrium. Due to Lp^E ^ the Liouville operator is not a self adjoint operator like in 
equilibrium. Therefore {nq\Ljq) ^ {Lnq\jq) and can therefore not simply be calculated. One 
possibility is to express the frequencies flnjit) and ^jj(t) with the help of a modified projec- 
tion operator p6[ by the one time structure factor, correlation functions and susceptibilities. 
The equation (|3^ ) can therefore be written in the form 

--^ = S,,{t) - Re[Sq{t) g2fi„,(0) + t)g2%(0) 

+iq dt'(f)Ut, t')Sr^{0)W^,{q; t' , t) (35) 

iqNjn{q;t',t)yMqit',t)] 





There appear 4 new memory functions Wmj and Umj with m = 1, 2 in this equation. 
Wmj, Umj are the correlation functions of the fluctuation of the variable ^4^(0) G {nq,jq) 
and the current of the fluctuating force LFq{t',t) between t' and t or the fluctuating force 
respectively, translated backwards in time to t = 0. 

Wmjit',t) = {ArnmeM-^Lt')LFqit',t)) (36) 

Umjit',t) = (A„(0)| exp(-zLt')F<?(t',t)) (37) 
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In equilibrium the memory function U vanishes exactly. The fluctuating force Fq{t',t) is 
perpendicular to Amit') and Umj(t',t) = {exp{+iLt')Am{0)\Fg{t',t)) = 0. The memoryfunc- 
tion Wnj is in equilibrium equal to the fluctuating force correlation function Mg. Together 
with the fluctuation dissipation theorem therefore the convolution integrals in time in ( |35D 
cancel. The flrst line in ( P5| ) vanishes, since in equilibrium Qjj = and q^^nj = Sjj/Sg. 
Equations (|35|j3^ , |26|) and (^) are an exact set of equations for one time and two time quan- 



tities, describing the change of the structure and the correlations and susceptibilities if e.g. 
by an initial quench a nonequilibrium probability density is imposed. 

III. APPROXIMATIONS 

To close the equations, closure relations for the memory functions Mg, Wmj, Umj and the 
correlation function of fluctuation forces Mg have to be found. There is no systematic way of 
calculating the memory function, since no small parameter is available for dense liquids close 
to the glass transition. The only guide available is the requirement, that the equations should 
reduce to the equ. ([Tl|) and (§) in equilibrium and that certain mathematical requirement 
intrinsic to the projection operator formalism are fulfllled. It is very plausible to express M 
as a quadratic form of the correlation functions since it has the structure of a correlation 
function of fluctuating forces. 

M{q;h,t) = ]-Sl{t)Y.v\q;k,q-k-t)(Pu{tut)(t)g-k{ti,t). (38) 

k 



The vertex f (g; k,q — k; t) is the same as in (12) if the direct correlation functions are 
understood as time dependent quantities. The memory function M describes the influence 
of a force fluctuation at t' on its surrounding integrated between t' and t. Since M cannot 
be written as a linear functional of the integrated response the next simple approximation 
possible is the following: 

M{q;ti,t) = Mg{t,t)- (39) 
^U'^^) E r dt'v\q- k,q-k- t'){-ik)N,n{k- 1', t)0g_fc(t', t) 
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where Njn{k; t', t) can be identified with the susceptibihty as explained above. In equihbrium 
this guarantees , that Mq = Mg. Therefore the two last integrals in eq. (|32| ) cancel for times 
ti = t. This property has to be fulfilled rigorously |^ . It can therefore be used to calculate 
M{t, t) selfconsistently. Although M and M are equal in equilibrium, the formal expressions 
in nonequilibrium are different. Choosing M = M has the consequence that a generalized 
fluctuation theorem of the second kind |^ would be fulfilled. It can be shown that this 



immediately excludes the possibility of the aging properties found e.g. in simulations |20 



^ . Only a transient aging phenomena caused by the time dependence of the static structure 



factors would be observed in this case. Another possibility to introduce approximations is 
first to integrate the third integral in ( ^2]) by parts. This leads to 

= -Mtut) (40) 
-S,(t)M(t,t)+0,(t,O)M(O,t) 

+ f'dt'(l)g{h,t')^it',t) (41) 
Jo 

+ f\-tqNj^{q;t',ti))Mq{t',t) 







with 



S(t',t) = Sjlit')^^^^ - ^qS-}it')Mgit',t) (42) 
With the approximation 

M{q;h,t) = Slih)J2v\q;k,q - k;t)Nj4k;h,t)(l){q - k;h,t) (43) 



the resulting structure of the theory is very similar to the one of the p - spin model [13 



This is also true for the theory introduced by the approximations (|38| , |39[) . Both reduce to 
EMCT in equilibrium. Which of the approximations is more appropriate for real liquids has 
to be tested and is part of future work. 

Similar approximations are possible for the new memory functions U and W. But since 
there are less mathematical boundary conditions to be fulfilled than for M and M there is 
more freedom in the choice of the approximation. Different approximation will lead only to 
quantitatively but not qualitatively different results for Sg. 
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IV. CONCLUSIONS 



With the help of non equihbrium projection operator techniques it is possible to derive 
an exact set of equations for one-time quantities like the time dependent equal time structure 
factor and two-time quantities like correlation functions and fluctuation propagators. The 
memory functions can be approximated as in the EMCT by the correlation functions and 
response functions. It is interesting to note that in the formulation presented here, a theory 
for the equal time structure factor (which is trivial in the p - spin model) and not the 
potential energy as in the p -spin model is derived. The structure factor is directly accessible 
to experiments and it would be especially interesting to study aging phenomena in colloidal 
hard sphere glasses as it is done for colloidal clays ||30|, where the structure factor can easily 
be measured within light scattering experiments. It is already obvious from the equations, 
that the structure factor lirrit^ooSqit) in the glass will not be the structure factor of the 
corresponding equilibrium liquid, but will depend on the nonergodicity parameters i.e. on 
\imt-,oo^iT^t^-*oo4>q{tw,t) via the memory functions Wjj,Unj and the correlation function 
of fluctuating forces M. The aging scenario to be expected from the equations presented 
here can be formulated within the bifurcation scenario of the glass transition derived from 
mode coupling theory 0. This picture can best be explained in discretized models, where 
only a discretized set of wave vectors is used. In this case the wavevector components of 
the equilibrium structure factor S^'^^ can be considered as control parameter of the mode 
coupling models. In the space formed by the components exists a critical hypersurface which 
separates the control parameter space in a liquid subset and a glassy subset. Since in the 
nonequilibrium case the value of Sg{t) for t —>■ oo depends on the nonergodicity parameters, 
the critical hupersurface will also depend on the two time correlation function. A quench can 
now be considered in the following way. By starting in the liquid subset the system, which 
can be considered as a a point in the control parameter space, starts to flow towards the 
critical hypersurface. As it approaches this surface relaxations become slower and slower. 
Due to the glass transition at the critical hypersurface the system can not reach the point, 
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which would correspond to the thermodynamic parameters in equihbrium. Instead it will 
very likely be trapped on the (nonequilibrium) critical hypersurface, which also depends on 
the nonergodicity parameters. In hard core like systems the glass transition is driven by 
the nearest neighbour repulsion. This can be inferred from the fact, that in the EMCT 
equations the glass transition occurs, if the first peak of the structure factor becomes big 
enough. Therefore trapping on the critical hypersurface will lead to a drastic change in the 
sensitivity of the peak height on the quench temperature (or density), once the system is 
quenched below Tc (above Uc). This qualitative aspect can be seen in simulations of binary 
Lennard Jones sytems p2[. In this simulation the area under the first peak of the radial 
distribution function g{r), which roughly corresponds to the height of the first peak of the 
structure factor Sg, does not seem to increase beyond a certain threshold, if the system is 
quenched below Tc. As a next step the details of the theory presented will be studied in 
simplified models before calculating the properties of realistic glass forming liquids. This 
will be part of future studies. 
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